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^ ' We consider the arrival time distribution defined through the quantum probability current for a Gaussian wave 
f"*'! \ packet representing free particles in quantum mechanics in order to explore the issue of the classical limit of 
arrival time. We formulate the classical analogue of the arrival time distribution for an ensemble of free particles 
represented by a phase space distribution function evolving under the classical Liouville's equation. The classical 
probability current so constructed matches with the quantum probability current in the limit of minimum uncertainty. 
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qq , Further, it is possible to show in general that smooth transitions from the quantum mechanical probability cur- 
rent and the mean arrival time to their respective classical values are obtained in the limit of large mass of the particles. 
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l.INTRODUCTION 



It is generally believed that a necessary requirement for the universal validity of quantum mechanics is that its 
results in the macroscopic limit must agree with those of classical mechanics, because the latter is well verified in the 
macroscopic domain. However, there exist vexed problems regarding the connection between classical and quantum 
mechanics; the question whether quantum mechanics in the macroscopic limit is completely equivalent to classical 
mechanics remains the focal point of diverging view points. This is poignantly reflected in the various controversies 
persisting in the relevant literature . Several naive interpretations of the classical limit of quantum mechanics 

based on approaches such as the h — > limit, the large quantum number (N — * oo) limit, or the Ehrcnfcst theorem, 
are all riddled with well known difficulties Einstein [2] and Pauli 3] strongly advocated the tenet that in the 

macroscopic limit, not only the localised wave functions but all physically admissible solutions of the Schrodinger 
equation must lead to predictions equivalent to those obtainable from classical mechanics. Such comparison between 
the two mechanics can be meaningful only within the framework of the ensemble interpretation. Thus the classical 
limit problem boils down to probing whether there is complete equivalence in the macroscopic limit between the 
empirical predictions of classical and quantum mechanics with respect to the properties of the same initial ensemble. 
This is the spirit which motivates the present investigation. 

For complete equivalence between classical mechanics and the macroscopic limit of quantum mechanics the fol- 
lowing criterion is necessary. In the classical limit all the measurable properties of a quantum mechanical ensemble 
corresponding to any normalizable wave function ty(x,t) should be equally reproduced by the classical phase space 
formalism using a distribution function D(x,p,t) utterly determined by the classical phase space description where 
the time-development of D(x,p,t) is in accordance with the classical Liouville's equation, and where the initial phase 
space distribution function for the ensemble of particles is taken matching with the initial quantum position and 
momentum distribution. In the current investigation we formulate the classical phase space distribution in a way 
which is completely classical unlike the one that is called the quantum phase space distribution function such as the 
Wigner distribution function [5| . The latter is essentially a quantum entity obtained by directly using the expression 
of the wave function, and is constructed to reproduce the results of quantum mechanics, but it does not satisfy the 
classical Liouville's equation. So, the Wigner distribution function, not being a positive definite quantity in general, 
does not provide the results of a classical phase space evolution. In contrast we formulate a phase space distribution 
function D(x,p,t) that is positive definite and also satisfies the classical Liouville's equation. The motivation for this 
work is to study the comparison between quantum mechanical results and those obtained from a purely classical phase 
space description by formulating a proper classical counterpart of the quantum ensemble. Here our focus is on the 
arrival time of the free particles but one can also investigate the quantum-classical comparison for other dynamical 
variables for particles in various types of potentials using the same approach. 

In recent years there has been an upsurge of interest in understanding the concept of time of arrival for a quantum 
particle Q|. In general, the issue of providing physically meaningful definitions of experimentally measured times in 
varied arenas such as tunneling times, decay times, dwell times, delay times, arrival times has gained importance jjjj]. 
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In this paper we are specifically concerned with the issue of arrival time. In classical mechanics, a particle follows 
a definite trajectory; hence the time at which a particle reaches a given location is a well defined concept. On the 
other hand, in standard quantum mechanics, the meaning of arrival time has remained rather obscure. Indeed, there 
exists an extensive literature on the treatment of arrival time distribution in quantum mechanics [8]. A consistent 
approach of formulating a definition for arrival time distribution is through the quantum probability current [9( . The 
quantum probability current if defined in an unambiguous manner contains the spin of a particle, as was pointed 
out by Holland 10(. Recently it has been shown using the explicit example of a Gaussian wave packet that the 
spin-dependence of the probability current leads to the spin-dependence of the mean arrival time for free particles 



111 ]. This effect, if experimentally observed, should place the probability current approach to mean arrival time on 
a firmer footing. A key issue for any definition of time of arrival in quantum mechanics is to secure an acceptable 
classical limit of the arrival time formulation. We formulate a classical analogue of the arrival time distribution for 
free particles obtained via the quantum probability current. Aspects of the quantum-to-classical transition for the 
arrival time distribution are then investigated. 

2. ARRIVAL TIME DISTRIBUTION 

We begin our analysis with the standard description of the flow of probability in quantum mechanics, which is 
governed by the continuity equation derived from the Schrodinger equation given by 

||*(x,i)| 2 +V.J(x,i)=0 (1) 

The quantity J(x, t)=^-(^V 1 &* — 4 r *V4 r ) defined as the probability current density corresponds to this flow of 
probability. We use this current to define the arrival time distribution for free particles. Interpreting the equation of 
continuity in terms of the flow of physical probability, the Born interpretation for the squared modulus of the wave 
function and its time derivative suggest that the mean arrival time of the particles reaching a detector located at 
x = X may be written as 

_ = / oo |J(x = X,t)|ftft 

T f Q °°\J(x = X,t)\dt U 
Henceforth, for simplicity we shall restrict ourselves to only one spatial dimension. One should keep in mind that 
the definition of the mean arrival time used in Eq( s).( 2) is not a uniquely derivable result within standard quantum 



mechanics. However, the Bohmian interpretation 12| of quantum mechanics in terms of the causal trajectories of 
individual particles implies the above expression for the mean arrival time in a unique and rigorous way [13I ]. It 
should also be noted that in ccratin situations J(A, t) can be negative over some time interval provided the initial flux 

n 

J(X,t = 0) is negative [14|. In order to account for the back flow effect in such cases, the decomposition of J(X,t) 
into right and left moving parts could be undertaken. However, our present analysis is carried out using a simple 
example that is free from such complications. 

The standard Schrodinger probability current defined through the continuity equation in non-relativistic quantum 
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mechanics, however, suffers from an inherent ambiguity since the continuity equation remains satisfied with the 
addition of any divergence free term to the current. This feature was exploited to formulate alternative causal 

n ... ..... n 

models |15| . Finkelstein has analysed the consequent ambiguities of arrival time distributions 16]. However, it 
was shown earlier by Gurtler and Hestenes that the problem of non-unique probability current doesn't exist for 
the relativistic Pauli theory for the electron if the probability current is inclusive of a spin-dependent term 1t| |. 
Holland demonstrated the uniqueness of the conserved probability current in the non-relativistic limit of the 
Dirac equation. This probability current differs from the standard Schrodinger probability current by the presence 
of a spin-dependent term which persists even in the non-relativistic limit. It has been further argued that the arrival 
time for a free particle computed using the unique probability current should exhibit an observable spin dependence 
111 ]. However, for the case of massive spin-0 particles it has been shown recently by taking the non-relativistic limit 
of Kemmer equation [psl ] that the unique probability current is given by the Schrodinger current, and hence, the 
Schrodinger current gives the unique probability current density or the unique arrival time distribution for spin-0 
particles [191 ]. In the present analysis we restrict our attention to massive spin-0 particles only. 



3. CLASSICAL-QUANTUM CORRESPONDENCE 

Let us now consider a Gaussian wave packet representing a quantum free particle moving in 1-D whose initial wave 
function 0) and its Fourier transform 3>(p, 0) are respectively given by 



*(a;,0) 



(2W) 1/4 x/TT~iC 



X ikX 4V(l + iC) } 



(3) 



2 2 y/4 rv^ (1+lC) | 



*<P.0)=(^j e I " 'J (4) 

where the group velocity of the wave packet u = hk/m = p/m. For generality we have taken the initial Gaussian 



wave function 0) which is not a minimum uncertainty state (AxAp — {h/2)yl + C 2 > H/2), but which could 

wave function ^(x, i), the quantum position probability 
density pQ(x,t) and the probability current density Jq(x,t) at a particular location x are then respectively given by 



*(x,t) 



(2W) 1 / 4 v /l + *(C+^) 



e *(x-,irt) exp . 



(x - utf 



4oo 2 



(5) 



p Q (x,t) = \^(x,t)\' 



exp 



(x - ut) 2 



2oV 



1 + (C+^) 2 



(6) 



Jc)(x,t) = pq{x,t) I u + 



2mco 2 



l + (^+») 2 



(7) 
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In order to elucidate the classical counterpart of the quantum probability current, we now construct a classical 
formulation of arrival time for an ensemble of free particles. We take the initial phase space distribution function 
for the ensemble of particles as a product of two Gaussian functions matching with the initial quantum position and 
momentum distributions from Eq(s).(3) and (4) as 

D ( Xo , Po ,0) = MxoM 2 l*( P o,0)| 2 = ^== e* P {- 2ao ^ +c2) - ^(Po-P)^ (8) 

where the variables xq and po are the initial positions and momenta of the particles. Note that our approach to 
compare the quantum and classical predictions is not contingent to any particular initial form of the wave function. 
The key point of this scheme is to choose the initial classical ensemble in such a way that it reproduces the initial 
quantum position and momentum distributions. Classically of course there are other choices for Dq(xq,p$, 0). But 
in quantum mechanics, due to the uncertainty principle, given a wave function ?p(x,t), the momentum space wave 
function <p(p, t) is automatically fixed by the Fourier transform of ip( x J t)- I n this wav t ne position probability density 
\ip\ 2 and the momentum probability density |</>| 2 are correlated in quantum mechanics. There is no such restriction for 
the position and momentum densities in classical statistical mechanics. But it is quite reasonable to take the initial 
classical phase space distribution exactly matching with the initial quantum position and momentum probability 
densities in order to compare the results obtained from the dynamical evolutions of classical and quantum mechanics. 
This is precisely the motivation to take the initial phase space distribution Do(xo 1 po, 0) in a way given by Eq(s).(8). 

Now to obtain the time evolved density function D{x 1 p 1 t) we focus on the classical dynamics of freely moving 
particles. The Hamiltonian is H = p 2 /2m and the Hamilton's equations are x = pt/m + xq and p = pa where the 
variables xq andpo are the initial position and momentum of the particle which are respectively given by xq — x—pt/m 
and po — p. Substituting these values of xq and po in the expression of Dq(xq,po, 0) we obtain the time evolved 
distribution function D(x,p, t). This is because here we are considering the free evolution of an ensemble of particles 
whose initial positions (xq) and momenta (pq) are distributed according to the initial density function Dq(xo 7 p$, 0). 
The time evolved phase space distribution is then given by 

D(x p t) I exp L (gzgg _ ^ 2 (p-p) 2 ] (9) 



At this stage it is instructive to write down the the Wigner distribution function which is calculated [20J from the 
time evolved wave function (ty(x,t)), and is given by 

1 f°° 

D w (x,p,t) = — V*(x + y,t)V(x-y,t) exp{2ipy/h}dy (10) 
By substituting the value of ^(x + y, t) and ^>(x — y, t) using Eq(s).(5) we obtain 
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FIG. 1: The position probability densities pg(x,t) and pc(x,t) are plotted for varying mass of the particles (in atomic mass 
units) with ao = 10~ 5 cm, u — 10 3 cm/sec, C=10 and t=10 -5 sec. 

Note here that the Wigner function Dw(x,p, t) is not identical with our classical phase space distribution D(x,p, t) 
where in the spirit of a completely classical description we have not included any position-momentum correlation. 

We consider a classical statistical ensemble of particles defined by the phase space density function D(x,p, t) in one 
dimension. Then the position and momentum distribution functions are respectively pc{x,t) — J D(x,p,t)dp and 
pc{p,t) — J D(x,p,t)dx. The classical position probability distribution for this ensemble is given by 



Pc(x,t) = J D(x,p,t)dp 



1 



(27r ( 7 2 ) 1 /2 v /l + C 2 + 



h 2 t 2 



exp ■ 



(x - utf 



2ao 2 (l + C2 + 3 ^) 



(12) 



4m 2 ao 4 - 

All the density functions are assumed to be normalized and D(x,p, t) satisfies the classical Liouville's equation 
given by 

dD(x,p,t) ,dD(x,p,t) ,dD(x,p,t) 



21] 



+ x- 



P- 



dt dx 

Since for free particles p = and x = p/m we have 

dD(x,p,t) p dD(x,p,t) 
dt m dx 

Integrating the above equation with respect to p one gets 

1 



dp 



= 



dpc(x,t) d_ 
dt dx 



p(x,t)p c {x,t) 







(13) 



(14) 



(15) 



where p — J pD{x,p,t)dp/ J D(x,p,t)dp is the ensemble average of the momentum. Defining v(x,t) — p(x,t)/m as 
the average velocity, one obtains 



dpc{x,t) 9 

^— + -j c (x,t) = 



(16) 



where Jc{x,t) and v(x,t) represent the mean motion of the continuum matter at (x,t). Eq(s).(16) is the equation of 
continuity for the continuous density function pc(x,t) of a statistical ensemble of particles. The expression for the 
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FIG. 2: The probability current densities Jg(x,t) and Jc(x,t) are plotted for varying mass of the particles (in atomic mass 
units) at a detector location X=10 cm with <to = 10~ 4 cm, u = 10 3 cm/sec, C=100. 



classical probability current density is given by 

Jc(x,t) 



1 

in 



pD(x,p, t)dp 



(17) 



and is related to the mean velocity by Jc{x,t) — pc(x,t)v(x,t). 

Now substituting the expression for the time evolved phase space distribution function D(x,p,t) from Eq(s).(9) 
in Eq(s).(17) we get the expression for the current density or the arrival time distribution at a particular detector 
location x=X for this classical ensemble of free particles given by 

(x — ut)h 2 t 



Jc(x,t) = pc{x,t) < u + 



(18) 



[h 2 t 2 + 4m 2 (To 4 (l + C 2 )] 

If we impose here the minimum uncertainty condition viz., C = then one can check from Eq(s).(6), (7), (12) and 
(18) that both pQ{x,t)=pc{x,i) and Jq(X,1) = Jc(X,t) hold, i.e., the classical and quantum probability currents 
are similar. Thus, if we take the initial phase space distribution function for the classical ensemble of particles as a 
product of two Gaussian functions matching with the initial quantum position and momentum distributions then the 
classical arrival time distribution exactly matches with the quantum one provided the minimum uncertainty relation is 
satisfied. But in general the quantum and classical distribution functions arc different when the minimum uncertainty 
condition is not satisfied (C ^ 0). 

Though Jq(X, i) and Jc(X,t) are in general not equal for C / 0, the large mass limits of both are the same. 
This is seen from Figures 1 and 2 where the probability distributions and the currents are plotted respectively for 
different masses. It is apparent that in the large mass limit quantum distributions reduce to the classical distributions. 
The mass dependendence in the arrival time distributions and also in the position probability densities (for both the 
quantum and classical case) arises from the spreading of the wave packet. 

We now compute the mean arrival time f by substituting the expressions for the quantum current in Eq(s) . (2) . [One 
should note that though the integral in the numerator of Eq(s) . (2) formally diverges logarithmically, several techniques 
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FIG. 3: The mean arrival time r is plotted against the mass of the particles (in atomic mass unit) at different detector locations 
X — 5.1 cm, X — 5.2 cm, and X — 5.3 cm. for C=10, ctq = 0.0001 cm, u — 10 cm/sec. 



have been employed in the literature 22| ensuring rapid fall off for the probability distributions asymptotically, so 
that convergent results are obtained for the integrated arrival time. Here we have employed a simple strategy of 
taking a cut-off (t = T) in the upper limit of the time integral with T = (X + 3<tt)/u where <jt is the width of the 
wave packet at time T. In other words, our computations of the arrival time are valid up to the 3<r level of spread 
in the wave function.] It is instructive to examine the variation of mean arrival time with the different parameters 
of the wave packet. In Figure 3 we have plotted the variation of f with mass at different detector locations, keeping 
the group velocity u and initial width <jq fixed. One sees that the mean arrival time calulated by using the quantum 
curCnt Jq(A, i) as the arrival time distribution asymptotically approaches the classical result in the limit of large 
mass. 

4. SUMMARY AND OUTLOOK 



To summarize, in this paper we have investigated the quantum-to-classical transition of the mean arrival time 
defined through the probability current. We have formulated the classical arrival time distribution from the phase 
space distribution for a classical ensemble of particles. The expression for classical probability current constructed 
by us matches exactly with the quantum probability current in the limit of minimum uncertainty. We note that the 
uncertainty condition is not a stringent requirement for the case of the initial classical distribution. Thus the classical 
arrival time distribution Jc(X,t) will in general be different from the quantum distribution Jq[X,t) if we do not 
impose the minimum uncertainty restriction on the initial distribution. This issue needs to be explored further in 
order to have a deeper understanding of the quantum-classical comparison of arrival time. However, in the present 
example that we have constructed, the quantum results for the probability current and through it the arrival time 
distribution, approaches the classical result in the large mass limit. A number of schemes 8, fJ 3 23 j have been 
suggested in the literature for calculating the arrival time distribution such as those based on axiomatic approaches, 
trajectory models of quantum mechanics, attempts to define and calculate the arrival time distribution using the 
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consistent histories approach, and attempts of constructing the time of arrival operator, etc. It might be worthwhile 
to investigate the quantum-classical correspondence of the arrival time distribution using these different approaches. 
Such studies, if undertaken extensively, are not only expected to throw light on the comparitive merits of different 
arrival time formulations, but could also be of relevance to the behaviour of mesoscopic systems where a great deal 
of experimental activity is presently underway 24 1 . 



Classically we know that a point particle with uniform motion will reach a particular location at a time which 
is independent of the mass of the particle and depends only on its uniform velocity. In the discussion of classical 
limit of quantum mechanics it is usually assumed that the peak mean position of the wave packet moves according 
to classical trajectory derived from the Ehrenfest theorem. It could be argued though that one should not expect to 
recover an individual classical trajectory when one takes the classical limit of quantum mechanics. Rather, one should 
expect the probability distributions of quantum mechanics to become equivalent to the probability distributions of an 
ensemble of classical trajectories. The current investigation is concerned about this particular approach to test the 
quantitative equivalence between the classical mechanical prediction and the prediction obtained in the macroscopic 
limit of quantum mechanics. Here, what we see is that the mean time of arrival of a freely moving quantum particle 
computed through the probability current depends on the mass of the particle even if its group velocity is fixed. So it 
turns out that the characteristic of mean time in this framework is different from that of mean position. The predicted 
mass dependence of mean arrival time is, in principle, amenable for experimental verification, and is a clear signature 
of the probability current approach to time in quantum mechanics [jj]. 
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